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ABSTRACT 



> 

^ We present our analysis of the long cadence Kepler data for the well-studied Blazhko 

star RR Lyr, gathered during the first two quarters of the satellite's observations and 
I/-) covering a total of 127 d. Besides being of great importance for our understanding of 

RR Lyrae stars in general, these RR Lyr data can be regarded as a case study for 
observations of bright stars with Kepler. Kepler can perform high-precision photometry 
on targets like RR Lyr, as the saturated flux is conserved to a very high degree. 
The Kepler data on RR Lyr are revolutionary in several respects. Even with long- 
cadence sampling (one measurement per 29.4 min), the unprecedented precision ( < 
mmag) of the Kepler photometry allows the study of the star's extreme light curve 
variations in detail. The multiplet structures at the main frequency and its harmonics, 
typical for Blazhko stars, are clearly detected up to the quintuplets. For the first time, 
photometric data of RR Lyr reveal the presence of half-integer frequencies, linked to 
a period doubling effect. This phenomenon may be connected to the still unexplained 
Blazhko modulation. Moreover, with three observed Blazhko cycles at our disposal, we 
observe that there is no exact repetition in the light curve changes from one modulation 
cycle to the next for RR Lyr. This may be due to additional periodicities in the star, 
or to transient or quasi-periodic changes. 

Key words: stars: oscillations - stars: variable: RR Lyrae - stars: individual: RR 
Lyr, KIC7198959 - techniques: photometric. 



1 INTRODUCTION 

The pulsating star RR Lyr has been studied for over a cen- 
tury now (Szeidl & Kollath 2000). In 1948 this field star 



became the eponym of a class of variables previously called 
"cluster variables" (see Smith 1995). RR Lyrae stars are low- 
mass, horizontal branch stars (post-red giant stage) that are 
burning helium in their core. They play a key role in as- 
trophysics for their use as standard candles, as examples 
of stars in an advanced evolutionary stage, and because of 
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their "simple" radial pulsations. RRLyrae stars have typi- 
cal periods of ~0.2d to ~ld, amplitudes in the optical of 
0.3 mag up to 2 mag, and spectral types of A2 to F6. They 
pulsate in the radial fundamental mode (RRab stars), the 
radial first overtone (RRc stars) and, in some cases, in both 
modes simultaneously (RRd stars). A few RR Lyrae stars 
are suspected to be pulsating in higher-order radial overtone 
modes. 

A large fraction of RR Lyrae stars shows the Blazhko 
effect (Blazhko 1907), a periodic amplitude and/or phase 
modulation with a period of typically tens to hundreds of 
times the pulsation period. The origin of the Blazhko effect 
remains a matter of controversy. Rotational modulation of 
the distorted radial mode invoked by a strong dipole-like 
magnetic field (Shibahashi 2000) appears to have been ruled 
out as a cause(Chadid et al. 2004; Kolenberg & Bagnulo 
2009). Several models proposing a resonance between the 
main radial mode and another mode have been proposed - 
with another radial mode (higher overtone: Borkowski 1980, 
Moskalik 1986), or with a nonradial mode of low spherical 
degree I (Van Hoolst et al. 1998; Dziembowski & Mizerski 
2004) - but their validity remains to be proven. Stothers 
(2006, 2010) proposes a scenario that attributes the Blazhko 
variation to variable turbulent convection caused by, e.g., 
transient magnetic fields in the star (that would be hard 
to detect). Thus far, all models presently proposed for the 
Blazhko effect have shortcomings in explaining the variety 
of features shown by modulated stars (Kovacs 2009). 

The RRab star RR Lyr is the brightest star, eponym 
and prototype of its class, and it shows very strong Blazhko 
modulation. It has been extensively studied over the course 
of the past century, through photometric (e.g., Szeidl & 
Kollath 2000), spectroscopic (e.g., Preston et al. 1965), and 
spectropolarimetric (e.g., Chadid et al. 2004) data. How- 
ever, several aspects of its pulsation remain poorly under- 
stood, and to accurately model the star, we have to take 
into account complex physics that we are only beginning to 
uncover. 

By a fortunate coincidence, RR Lyr (a( J2000): 19 h 25 m 
28 s and S( J2000): +42° 47' 04") lies in the field of the Kepler 
space telescope (Kepler ID: KIC7198959), though its bright- 
ness is well above the saturation limit for the Kepler CCDs. 
The unprecedented quality of the Kepler data of the star 
leads us to several new findings that spur further modelling. 
To explore theoretical models of RR Lyr, it is necessary to 
know the physical parameters of the star with the highest 
accuracy possible. For this reason, a self-consistent atmo- 
spheric and abundance analysis from high-resolution, high 
signal-to-noise spectra of the star, with the aim of param- 
eter determination, was recently published by Kolenberg et 
al. (2010b). 

According to the General Catalogue of Variable Stars 
(Samus 2004) its V brightness changes within the range 7.06 
- 8.12 mag and its period is P = 0.56686776 d or about 13 
h 36 min, corresponding to a frequency of fo — 1.7640799 
d _1 . Kolenberg et al. (2006) found a (mean) frequency of 
fo = 1.764170 ± 0.000005 d" 1 from data gathered in 2003 
and 2004. 

In Section 2 we present the new observations of RR Lyr 
obtained with the Kepler satellite and describe the chal- 
lenges of the data reduction. Section 3 presents the data 
analysis and its results. In Section 4 we discuss our results 



in the framework of the current understanding of RR Lyrae 
pulsations and the Blazhko effect. Finally, some concluding 
remarks and the outlook for future investigations are given 
in Section 5. 



2 OBSERVATIONS 

The Kepler space telescope (Borucki et al. 2010) was de- 
signed to detect transits of Earth-like planets around Sun- 
like stars. Technical details on the mission can be found in 
Koch et al. (2010) and Jenkins et al. (2010a, b). In order 
to exploit the potential of the Kepler asteroseismic data 
(Gilliland et al. 2010), the Kepler Asteroseismic Science 
Consortium (KASC) was set up, grouping more than 350 
scientists from all over the world. KASC working group 13 
is dedicated to the study of the Kepler RR Lyrae stars. At 
the time of writing, over 30 RR Lyrae stars have been iden- 
tified in the Kepler field and are now observed by the satel- 
lite. About half of the Kepler RR Lyrae targets have turned 
out to be modulated (Kolenberg et al. 2010a; Benko et al. 
2010). With Kepler's unprecedented precision, we will find 
additional clues to solving the Blazhko riddle. 

The Kepler Mission offers two options for observa- 
tions: either long cadence (29.4 minutes) or short cadence 
(1 minute). The spacecraft "rolls" every three months to 
allow for continuous illumination of Kepler's solar arrays. 
RR Lyr was observed in long-cadence during the first roll 
of the Kepler survey phase between HJD 2454964.0109 and 
HJD 2454997.4812 (2009 May 12 - June 14; Ql data, 33.5 d, 
1628 useful data points), and during the second roll be- 
tween HJD 55002.01748 and HJD 55090.96492 (2009 June 19 
- September 16; Q2 data 88.8 d, 4097 useful data points). 
The total time span of the data set is nearly 127 d, more than 

3 complete Blazhko cycles, as can be seen in Fig.[l] There 
are a few small gaps in the data set, due to unplanned safe 
mode, loss of fine pointing events as well as regular data 
downlink periods. Note that except for around those gaps, 
nearly every single Kepler measurement is of excellent qual- 
ity. Long-cadence sampling is just sufficient to obtain a good 
coverage of the light curve, but rapid changes and "glitches" 
in the light curve are missed by this sampling, and can be in- 
vestigated with forthcoming short-cadence Kepler data. As 
a consequence of the sampling rate, the Nyquist frequency 
lies at 24.5 d" 1 . 

Fig. [I] shows the quality of the Kepler data on RR Lyr 
compared to one of the most precise published data sets on 
the star obtained from ground-based observatories (Kolen- 
berg et al. 2006). Despite the effort of organizing a multi-site 
campaign and combining all the standardized observations, 
the latter data set had a point-to-point scatter of over 0.005 
mag at best, and was characterized by nightly and weather- 
related gaps, typical for Earth-based observations. As a con- 
sequence, the Fourier spectra of such data sets obtained from 
the ground have a noise level 10-50 times higher than Kepler 
Fourier spectra and are subject to aliasing. 

2.1 Early Kepler photometry of RR Lyr 

Kepler's CCDs saturate between Kepler magnitudes 11 and 
12, depending on which CCD a target falls, as well as the 
target location on that CCD. This saturation spills up and 
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Figure 1. Comparison of the ground- based RR Lyr data, gathered from six different observatories, published by Kolenberg et al. (2006) 
and the Kepler Ql and Q2 data of the star transformed to the magnitude scale (top panels). Individual cycles are shown in Fig. [3] Bottom 
panels: Fourier transform of the data; the insert shows the window function. 



down the CCD column, and saturated flux is conserved to a 
very high degree (see Fig.[2|. Therefore, Kepler can perform 
high-precision photometry on saturated targets like RR Lyr. 
Due to bandwidth constraints, Kepler downlinks pixels se- 
lected for each target star (Haas et al. 2010). These pixels 
are selected algorithmically to maximize the signal-to-noise 
ratio for the target (Bryson et al. 2010), based on stellar 
properties in the Kepler Input Catalog (KIC) and various 
models of sky and spacecraft characteristics, including satu- 
ration behavior. In the case of RR Lyr, the Kepler magnitude 
in the KIC did not reflect maximum brightness. This, com- 
bined with the early and approximate state of the saturation 
model at the time of Ql and Q2 observations, resulted in the 
size of RR Lyr's saturation being underestimated. Therefore, 
not enough pixels were allocated to RR Lyr to fully capture 
all flux at its brightest maxima: the central saturation col- 
umn exceeded the assigned aperture. In Ql the consequent 
loss of flux was relatively small but in Q2 the loss of flux at 
maxima was significant. 

Fortunately, columns adjacent to the central saturation 
column also saturated, though to a lesser extent, and these 
columns were fully captured in all of Ql and Q2. This al- 
lows us to use those measurements where the central column 
is fully captured (away from maxima) to determine the ra- 
tio of the central columns flux to the flux in the adjacent 
columns, and use this ratio as a predictor of the flux in the 
central column when the central column is not fully cap- 
tured. This ratio, ro, is sensitive to various known pointing 
and focus systematics, so we determine ro for every measure- 
ment, and detrend ro using a robust piecewise-polynomial 
fit to those measurements where the central column's flux is 
fully captured. The resulting polynomial trend rt was con- 
sistent with known motion and focus systematics. We define 



a correction c — r t /ro which multiplies the central column's 
flux for those data points in which some of the flux falls out- 
side RR Lyrs pixels, as determined by the detrended ratio 
rd falling below a cutoff value. In Ql only minor corrections 
on the order of 10 per cent were required for the central 
column flux, and these were at the beginning and end of Ql 
(from MJD 54964.25 to 54972.76 and 54995.47 to 54997.17). 
In Q2 corrections to the central column's flux were required 
for every maximum and ranged from 30 to 60 per cent. The 
uncertainty in the correction is estimated, from the variance 
of rd, to be 0.45 per cent in Ql and 0.36 per cent in Q2 
(neglecting six measurements with large outliers in Q2). 

Once the central column flux is corrected, it is summed 
with the other columns in the aperture to create the esti- 
mated flux correction for each measurement. This flux cor- 
rection is applied to the output of the Kepler Science Op- 
erations Center pipeline, which corrects light curves by re- 
moving cosmic ray hits, background and correlations with 
known systematics (Jenkins et al. 2010a,b). 

An example of the original and corrected light curves 
in Q2 is shown in Fig. [3] In Ql the final flux corrections 
were on the order of at most 5 per cent. In Q2 the flux cor- 
rections range from 15 to 35 per cent. The uncertainty in 
the total corrected flux is about 0.25 per cent in Ql and Q2 
for measurements in which the correction was applied, com- 
pared with the flux uncertainty of 8 x 10 -6 in measurements 
where the correction was not applied because the flux was 
completely captured. 

2.2 RR Lyr's Kepler Magnitude 

Kepler processing does not provide calibrated Kepler mag- 
nitudes Kp. We produced a rough estimate of RR Lyr's 
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Figure 2. RR Lyr's downlinked pixels during quarter 2 (Q2), 
showing the saturation spill in the three central columns. The first 
cadence in Q2 is shown, when RR Lyr was near a minimum. For 
brighter cadences the central column saturation exceeded these 
pixels, but the two adjacent columns remained completely cap- 
tured. Note the differing row and column scales. 



ture and ambiguity is only in tying the RR Lyr light curve 
to the absolute Kepler photometry. 

As can be seen in Fig.[l]s upper right panel, RR Lyr's 
total light curve amplitude varies significantly over the 
Blazhko cycle. The total variation (from lowest minimum 
to highest maximum) spans 0.86 mag in Kp, with 0.84 mag 
for the largest amplitude light curves, and 0.47 mag for the 
smallest amplitude light curves. Hence, the total reduction 
of the light curve amplitude due to the Blazhko effect is 
about 44 per cent as determined from our Kepler data of 
RR Lyr so far. These values will be refined with the short 
cadence data to come. 

The GCVS lists a 7.06 - 8.12 (1.06 mag; see also total 
amplitude in Fig.[l]s left upper panel) variation in V magni- 
tude. Hence, with the wide (white-light) passband of Kepler 
we observe an amplitude reduction with respect to the V 
magnitudes by a factor of 1.23. The same factor is found 
when comparing the amplitude Ai of fo (the first Fourier 
component) in the V data published by Kolenberg et al. 
(2006) to the one obtained from our data (Table[l]). 

The possible deviations from the "real" Kepler mag- 
nitudes introduced by the corrections and their associated 
uncertainties have a negligible effect on the frequency values. 
For Q2, the quarter covering more than two Blazhko cycles, 
and for which the largest corrections were needed, the fre- 
quency values derived from the non-corrected and corrected 
data are identical to within the error bars. For the Ql data 
this is also the case but the comparison is less convincing 
since they cover less than one Blazhko cycle. 




MJD 

Figure 3. Comparison between the captured flux from RR Lyr 
(dashed line) and the corrected flux (full line) during Q2, with 
±1 cr error bars. 



magnitude in Ql and Q2 using 5 nearby quiet (95 per cent 
variation under 1 part per thousand) stars of known Kp. 
Because this method did not account for colour or uncer- 
tainties in the comparison star's magnitudes and measured 
flux, the Kp magnitude values here are only approximate. 
We emphasize, however, that the relative precision within a 
given quarter is unprecedentedly high, the approximate na- 



3 DATA ANALYSIS AND RESULTS 

A detailed visual inspection of the Kepler data of RR Lyr 
has already revealed changes between consecutive Blazhko 
cycles, as well as alternating higher and lower light curve 
maxima at certain phases in the Blazhko cycle (see fig. 4 
in Kolenberg et al. 2010a). Because of this, we can expect 
a more complex frequency spectrum than the "classical" 
"main frequency + harmonics + equidistant multiplet struc- 
tures" predicted for a monoperiodic star undergoing ampli- 
tude and phase modulation (Szeidl & Jurcsik 2009; Benko et 
al. 2009). Moreover, the frequency solution we obtain partly 
depends on the coverage, as long-term modulation, secular 
trends, or transient phenomena cannot be fully captured by 
the Fourier techniques we adopt. The Fourier analyses pre- 
sented here were performed with Period04 (Lenz & Breger 
2005). 

Without adjustments, there is a small shift between the 
Ql and Q2 data even after the calibrations described in Sec- 
tion[2] To create the best match between the two data sets 
(as shown in Fig.[T|, we fitted the complete frequency solu- 
tion (described in the next sections) to each set separately 
and subsequently shifted by the zero point of the obtained 
fit. We used the combined data set (Ql and Q2) to deter- 
mine frequencies, as a longer time base and better coverage 
of the Blazhko cycle yields more accurate (mean) frequency 
values. Though we do not have sufficient information to re- 
liably combine the Ql and Q2 data regarding their ampli- 
tudes, the obtained frequency values are not significantly 
influenced by the small zero point shifts. Subsequently, we 
applied the harmonic fit with the frequencies derived from 
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Ql and Q2 to the longest data set (Q2). The results are sum- 
marized in Table [T] and discussed below. Optimum values 
for the frequencies, amplitudes and phases were obtained by 
minimizing the residuals of the fit. Errors were determined 
through Monte Carlo simulations. 



3.1 Typical Blazhko patterns 

RR Lyr's high-amplitude (nonlinear) pulsation gives rise to 
a non-sinusoidal light curve, described in Fourier space by 
contributions from the main frequency fo and its harmon- 
ics. The Blazhko modulation manifests itself as equidistantly 
spaced multiplets around the main frequency and its har- 
monics (Szeidl & Jurcsik 2009; Benko et al. 2009). All these 
frequencies were found through successive prewhitening, and 
the triplet components are strongly present in the data. 
Therefore, in a first stage, we fit the data with a Fourier 
sum of the form: 



f{t) = A + J2l =1 [A k sin(27r(fe/ t + fa)) 
+A+ sin(27r((fc/o + h)t + </>+)) 
+A-sin(27r((fc/o-/ B )t + <0] 
+Bosin(27r(/ B t + (?f>s)), 



(1) 



From the combined Qf and Q2 data, we find /o=l. 76416 ± 
0.00001 d" 1 for the main frequency, and its detected har- 
monics are significant up to the 14th order. From the side 
peaks (triplet components) we find /b=0.0256 ± 0.0002 d _1 
for the Blazhko frequency, or a Blazhko period of 39.1 ± 0.3 
d. 

From the Q2 data alone we find / = 1.76422 ± 0.00001 
d"\ and / B = 0.0252 ± 0.0001 d" 1 , or a Blazhko period of 
39.6 ± 0.3 d. As mentioned above, we also fit the Q2 data 
(Table [T| with the values derived from Ql and Q2 because 
of the longer time base. 

The main pulsation frequency we obtain is an average 
value of the actual, varying pulsation frequency over the 
time interval considered. In fact, the amplitude and period 
of RR Lyr's pulsation change significantly over the star's 
Blazhko cycle(s), as is shown in Fig.[4] obtained with the 
analytical signal method (Kollath et al. 2002). The pulsation 
period varies between 0.5652 d (1.7693 d _1 ) and 0.5699 d 
(1.7547 d" 1 ). This corresponds to a period change (5P/P) 
of about 0.83 per cent, or nearly 12 minutes. The period 
change is not identical from one Blazhko cycle to the next 
(see thick line in Fig. [4] for consecutive Blazhko cycles). The 
Blazhko frequency may in fact also be variable, and the value 
we obtain depends on the coverage (see Section 4.1 1. Note 



that the harmonics of the Blazhko frequency (2/b, 3/b, etc.) 
do not reach the significance level in this data set, as they 
do in several other Blazhko stars (see, e.g., Chadid et al. 
2010). This may change with additional data of the star. 



3.2 Quintuplet components 

After subtracting a fit of the form of Equation]!] we see many 
additional peaks in the spectrum. First of all, we see clear 
evidence for quintuplet frequencies. Fig. [5] shows the resid- 
ual spectrum after prewhitening with the main frequency, 
its harmonics and the triplet components (top panel), and 
zooms around some of the quintuplet components (bottom 
panels). The quintuplet components are added to Tablejl] 
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Figure 5. Residual spectrum of the Kepler RR Lyr data af- 
ter prewhitening with the main frequency, its harmonics and the 
triplet components (top panel). Zooms around some of the quin- 
tuplet components are shown in the bottom panels. 



They become more easily discernible in the frequency spec- 
trum with increasing order. At low order, e.g., around fo 
and 2/o, other close frequencies have higher amplitudes (see 
Fig.[5|. There are no clear indications for the occurrence of 
multiplet components of order higher than the quintuplet in 
the Fourier spectrum of the Kepler data of RR Lyr so far. 



3.3 Period doubling 

It is not only around the main frequency and its harmonic 
components where we see additional peaks. After prewhiten- 
ing with the main frequency fo = 1.76416 d - , its harmonics 
and the significant multiplet components, we clearly detect 
other significant frequencies in the Fourier spectrum. Fig.[6] 
shows a zoom in the frequency range [1.5 - 10] d _1 (fo ~ 
5/o)- We clearly see additional signal around the half-integer 
multiples of fo- The occurrence of such half-integer compo- 
nents is connected to the period doubling phenomenon (see 
Section 4.2 l. As a consequence of period doubling we see the 



previously mentioned alternating higher and lower maxima, 
as can be seen in Fig.[4](most striking in the middle panel), 
in fig. 4 of Kolenberg et al. (2010a) and in figs 1, 2 and 4 of 
Szabo et al. (2010). 

The highest peak is detected near 3/2/o, followed by 
5/2/o and l/2/ . The frequency at 2.6643 ± 0.0026 d _1 is 
the highest among the peaks around 3/2/o. Its ratio to the 
main frequency is 1.5 10 ±0.001, so there is a significant devi- 
ation from the exact half-integer ratio. We also observe sev- 
eral other frequency peaks around the 3/2/o component, and 
around the other half-integer frequency components, very 
close to each other, as shown in Fig.[6] They form asymmet- 
ric bunches that are more peaked towards the right (higher 
frequency) side. 
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Figure 4. Relationship between the period change and amplitude modulation over the Blazhko cycle for RR Lyr. The crosses show 
the Kepler data for the three monitored Blazhko cycles, the connecting line our multi-frequency fit to the data. The thick line shows 
an instantaneous period determined using the analytical signal method (units to the right side of the panel). The period is longest just 
before the lowest amplitude phases and shortest just after the highest amplitude phases. Note the small changes in the period variation 
between consecutive Blazhko cycles, and the alternating high and low maxima at certain Blazhko phases. 
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Figure 6. Fourier spectrum of the Kepler RR Lyr data in the range [1.5 — 10] d , after subtraction of the main frequency, its 
harmonics, and the triplet and quintuplets components. The positions of the mean pulsation frequency and its harmonics kfo, and that 
of the half-integer frequencies 2fc ~ 1 /o, are indicated by vertical lines. 



3.4 Additional peaks 

After prewhitening with the main frequency, its harmon- 
ics, the observed triplet/quintuplet side peaks and the half- 
integer components, we observe many residual peaks around 
their positions (see Fig.[6| . Their spacing to the prewhitened 
frequency peaks is not equal or suspiciously close to the 
Blazhko frequency, so we cannot directly connect them with 
the Blazhko modulation. A changing Blazhko effect, due to 
longer cycles, secular trends, or transient phenomena would 
give rise to frequencies close to the typical pattern. Further 
Kepler data to come will clarify what is the origin of the 
additional peaks. 



3.5 Light curve modulation 

3.5.1 Modulation components 

The properties of the modulation components, specifically 
the triplet components kfo + /b and kfo — /b, in the fre- 
quency spectra of a Blazhko star can provide constraints 
for the models proposed to explain the modulation. Al- 
cock et al. (2003) found that the relative amplitudes of 
the first order modulation components are usually in the 
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Figure 7. Amplitude variation versus phase of maximum light. 
The big dot marks the beginning of the trajectory. 

range 0.1 < Af/A± < 0.3. From the Kepler data we find 
At/Ai = 0.267 ± 0.010 and A~ /A 1 = 0.049 ± 0.007. 

Table[2]lists the amplitude ratios, phase differences and 
their errors for the first ten modulation component pairs, 
defined as Rk — A^ /A^ and A(j>k = <p^ — 4>^ , with A^ and 
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Table 1. Amplitudes and phases of the pulsation and modulation frequency components of RR Lyr for the best fit to the Q2 data. The 
values displayed in italics correspond to combination frequencies not exceeding a signal-to-noise level of 3.5. The residuals of the fit to 
the data are 0.0127 mag. We also list the most prominent additional period doubling components that arc found in the data set. The 
uncertainty on /q is 1 X 10~ 5 d _1 and the uncertainty on fg is 2 X 10~ 4 d _1 . 
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<jy^ defined by Eq. 1. We also list the so-called asymmetry 

parameter Q k = — \ — -It , introduced by Alcock et al. (2003) 

A k +A k 

to quantify the degree of asymmetry in the peaks. The dis- 
tribution of the Q parameter for the Blazhko stars from the 
MACHO data base (Alcock et al. 2003) peaks at +0.3. The 
positive Q points to an asymmetry with higher amplitudes 
at the higher frequency lobes in the triplets, as is mostly 



the case in Blazhko stars. Finally, we list the power differ- 
ence of the side peak amplitudes, AA k — (A^) 2 — (A^) 2 . 
As Szeidl & Jurcsik (2009) showed in their mathematical 
description of periodically modulated sinusoidal oscillation, 
this quantity depends on the phase difference between am- 
plitude and phase modulation components. Hence, it is the 
more physically meaningful quantity to measure the asym- 
metry of the triplet than the amplitude ratios Rk- If the 
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Figure 8. Left panel: amplitude ratios A^/Ai, Aj_ /A^ , and A k /A-^ of the detected harmonic and triplet frequencies in RR Lyr, with 
error bars. Right panel: amplitude ratios A^~*~ /A~^~*~ , and A^~ /A^~~ of the detected quintuplet frequencies in RR Lyr, with error bars. 




-0.6 
-0.4 

-0.2 

0.2 
11.4 







4 "(1=0-679 


A 

/ \ 






1 







-0.6 
-0.4 

-02 

02 
0.4 





1/1=0.085 










!f=O.I70 















ifc=0.509 
















k 


V=0.255 

A" 




Ma. 








\^ 






1 0.5 


1 1.5 




|(=0.594 









i «M>.934 


if 









Figure 9. Variation of RR Lyr's Kepler light curve (magnitude versus phase) for the second full Blazhko cycle in Q2, showing 12 
snapshots (connected points) roughly equally spaced in time; each sixth pulsation cycle is plotted. Blazhko phase tp = corresponds to 
maximum light amplitude. The dots are all the data in the Blazhko cycle folded with the main period. 



phase difference between the amplitude and phase modu- 
lations lies between and 7r, then > and the plot 
showing amplitude variation versus phase of maximum light 
has an anticlockwise progression, and vice versa (see Fig.[7|. 
Table|] shows, through its values of Rk > 1 and AA|, that 
the triplet components in RR Lyr always show asymmetry 
to the right side. However, the asymmetry of the side peak 
amplitudes decreases towards higher orders, reflected in in- 
creasingly smaller Rk and AA\ values. Also their phase dif- 
ference A(f>k, after peaking at order k = 2, decreases towards 
higher orders. 

Jurcsik et al. (2005) pointed out that generally the side 
peaks show a less steep decrease in amplitude towards higher 



orders than the successive harmonics of /q. Fig.|8](left panel) 
shows the amplitude ratios A^/Ai, A~^ /A^ , and A^ /A~^ . We 
confirm the less steep, though not "linear" decrease for the 
modulation components in RR Lyr. The decrease at higher 
orders is less steep for the left sidepeak components. For the 
first time, we also show a similar figure for the quintuplet 
components (Fig,|8j right panel). Their decrease, especially 
at higher orders, is similar to that of the triplet components. 



3.5.2 Fourier parameters 

RR Lyr's light curve dramatically changes its shape over the 
Blazhko cycle. Fig. [9] shows the variation of RR Lyr's light 
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Table 2. Side lobe amplitude ratios R^, phase differences A<f>k (x2tt), asymmetry parameters Qk, and power differences AA| as defined 
in the text, and their respective errors (based on Monte Carlo errors) for the Kp data of RR Lyr. k denotes the multiplet order. Values 
with a significance below 3ct are given in italics. 
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curve, through 12 snapshots evenly taken over the Blazhko 
cycle. 

We applied a time-dependent Fourier analysis (Kovacs 
et al. 1987) to derive the Fourier parameters for subsets of 
data. The data set was subdivided into segments consisting 
of either one pulsation cycle (26-28 points in a group), or 
four pulsation cycles (109-111 points in a group), where the 
pulsation cycle was defined starting from minimum bright- 
ness. The data were phased with a constant period, the av- 
erage period P — 0.566843 d obtained from our frequency 
analysis. A seventh-order harmonic fit, applied to these sub- 
sets of data, proved to yield the best fits for our sampling. 
The so derived Fourier parameters, i.e., the amplitude ratios 
Rki = Ai c /A 1 and the epoch-independent phase differences 
fki = fk^kipi, offer a way to quantify the shape of the light 
curves. We also used the analytical signal method (Kollath 
et al. 2002) to analyze the data, yielding consistent results, 
but due to the gaps in the data we used the time-dependent 
Fourier analysis to present the continuous variation over the 
whole Blazhko cycle. 

Time-dependent Fourier parameters derived from seg- 
ments four periods long show a smooth time variation re- 
vealing the Blazhko time scale; this procedure averages out 
the period-to-period oscillation caused by the period dou- 
bling phenomenon which is seen when segmenting on single 



periods (Fig. 10 1. 

To show the light curve changes over the Blazhko cy- 



cle we use the four-cycle groups (Fig. Ill, thus averaging 
out the period doubling effect. The interrelations between 
Ai and R21 (amplitude - amplitude ratio) and Ai and 9521 
(amplitude - phase difference) describe loops with particu- 
lar shapes, which models for the Blazhko effect should be 
able to reproduce. The A\ - R21 interrelation shows a dou- 
ble loop, of which the largest part has a counter-clockwise 
progression. For Ai - <p2i the (single) loop is wider towards 
lower Ai values, and its progression is clockwise. The loops 
for consecutive Blazhko cycles do not completely overlap: 
there is no exact repetition from one Blazhko cycle to the 
next. 

Fig. [12] shows the variation of the Fourier parameters 
during the modulation. The variability of A\ is almost si- 
nusoidal, which shows that the modulation itself is not very 
nonlinear. R21 reaches its maximum when Ai is at minimum. 
This is an observation that imposes constraints on the the- 
oretical models for the Blazhko effect. Generally, stars with 
smaller amplitudes have more sinusoidal light curves. If dur- 
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Figure 10. Comparison of the Fourier parameters (FP) Ai de- 
rived from one-cycle and from four-cycle segments. The effect of 
period doubling is clear in the parameters derived from one-cycle 
segments. The oscillation of the parameter A\ reflects the alter- 
nating height of successive light curves. When taking four-cycle 
segments this effect is smoothed out. 



ing the modulation pulsation energy is removed from the 
mode, the star is expected to behave as if being excited to a 
lower amplitude: its variability should be more sinusoidal - 
closer to the linear regime - and R21 should decrease. We ob- 
serve the opposite: the pulsation is more nonlinear at smaller 
amplitudes. This may suggest that the modulation is caused 
by a mode resonance. Detailed hydrodynamical models are 
needed to investigate the origin of this nonlinearity. 



4 DISCUSSION 

4.1 Blazhko modulation 

For many decades, RR Lyr's Blazhko period was known to 
be 40.8 d (see, e.g., Szeidl & Kollath 2000). Kolenberg et 
al. (2006, their table 6) list a summary of the Blazhko peri- 
ods of RR Lyr quoted in previous papers devoted to the star. 
Unfortunately, error bars are rarely given and hard to recon- 
struct. In the course of the past decades, the Blazhko period 
appears to have shortened. From the Ql and Q2 Kepler data 
we find a value for the Blazhko period Pb = 39.1 ± 0.3 d, 
corresponding to a Blazhko frequency of /b = 0.0256 d _1 . 
From the Q2 Kepler data alone, we find Pb = 39.6 ± 0.3 d, 
corresponding to a Blazhko frequency of /b = 0.0252 d _1 . 
It is well established that periodic amplitude and phase 
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As shown by Jurcsik et al. (2009), the mean parame- 
ters of a Blazhko RR Lyrae star vary over the course of the 
Blazhko cycle, such as the star's mean temperature, luminos- 
ity and radius. Due to these changes, also the star's intrinsic 
period, determined from the Kp light variations, varies be- 
tween 0.5652 d and 0.5699 d (Fig.|2| for the Q2 data. As the 
Blazhko modulation appears to vary from one Blazhko cycle 
to the next, it is likely that the Blazhko period itself is also 
variable. 
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Figure 12. Variation of the Fourier parameters during the 
Blazhko modulation. 



modulation results in equidistantly spaced multiplets (Szeidl 
& Jurcsik 2009; Benko et al. 2009). Though RR Lyr is one of 
the best studied stars of its class, it is the first time a signif- 
icant detection of quintuplet components is reported in the 
star. Since the first observation of quintuplet components 
by Hurta et al. (2008) in RV UMa, higher-order multiplet 
components (quintuplets and higher) have been found by 
several authors in high-quality data of Blazhko stars, such 
as MW Lyr (Jurcsik et al. 2008), SS For (Kolenberg et al. 
2009) and the CofloTtarget V1127 Aql (Chadid et al. 2010). 



4.2 Period doubling 

With the Kepler photometry, the period doubling phe- 
nomenon in RR Lyrae stars was reported for the first time 
(Kolenberg et al. 2010a). Even for a star as well-studied as 
RR Lyr, it had never been detected before. This is not only 
due to the much higher accuracy of the Kepler data, but 
also to the fact that from the mostly single-site observations 
of the star, consecutive pulsation cycles cannot be observed 
for stars with periods of typically half a day. The alternating 
higher and lower maxima would be easily missed by single- 
site observations due to daily gaps. Also recent multi-site 
observing campaigns of the star (e.g., Kolenberg et al. 2006) 
did not lead to a detection of the effect. 

The period doubling and alternating heights of max- 
ima are securely detected in two other Kepler stars, namely 
V808 Cyg (KIC4484128) and V355 Lyr (KIC755345) (see 
also Szabo et al. 2010). A further four Blazhko RR Lyrae 
stars in the Kepler field possibly exhibit the effect: V2178 
Cygni (KIC3864443), V354 Lyrae (KIC6183128), V445 
Lyrae (KIC6186029) and V360 Lyrae (KIC9697825) (see 
Benko et al. 2010). So far the period doubling phenomenon 
has only been found in Blazhko stars and, though not all 
of our Kepler Blazhko stars seem to show it, there may be 
a connection between its occurrence and the Blazhko effect 
(Szabo et al. 2010). The period doubling does not maintain 
the same strength along the Blazhko cycle. In some parts 
of the Blazhko cycle, the alternating light curve shapes are 
very obvious, in others indiscernible (see Fig.[4|. 

It is also of interest that period doubling occurs in 
RV Tauri stars and in many hydrodynamical models of Pop. I 
(e.g., Moskalik & Buchler 1991; Aikawa 2001) and Pop.II 
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cepheids (e.g., Buchler & Kovacs 1987; Buchler & Moskalik 
1992). 

We refer to Szabo et al. (2010) for a more detailed study 
of the period doubling phenomenon from the Kepler data. 
They also investigate and simulate the deviation of the ratio 
from exact half-integer values. The latter is due to a com- 
bined effect of the Blazhko modulation and the temporal 
onset and disappearance of the half-integer frequencies They 
suggest that a 9:2 resonance between the radial mode and 
the 9th-order radial overtone is responsible for the period 
doubling. 

Moskalik & Buchler (1990) studied the period doubling 
found in Pop. I and Pop. II cepheid models in detail. Through 
the analysis of the limit cycle stability they showed that the 
3:2 (in Pop. I) and 5:2 (in Pop. II) resonances are responsible. 
They demonstrated the coincidence of the period doubling 
and the resonance center, and provided a description of the 
mechanism in which the resonance destabilizes the limit cy- 
cle and causes the bifurcation. They showed that every half 
integer resonance, in principle, is capable of causing period 
doubling. The detailed modeling of the period doubling in 
Blazhko RR Lyrae stars will be discussed in a forthcoming 
paper by Kollath et al.. 

4.3 Residual frequencies and additional cycles 

The Kepler data of RR Lyr indicate that there is no strict 
repetition from one Blazhko cycle to the next, as is also 
reflected in Figs. 5 and |11| The meticulous efforts to recon- 
struct the light variation of RR Lyr described in Section |2~2| 
corroborate the reality of this variation. In their recent high- 
quality CoRoT data of VI 127 Aql, Chadid et al. (2010) 
found a small but significant shift of the maximum phase 
over five consecutive Blazhko cycles. On a longer time scale, 
irregularities in RR Lyr's Blazhko cycle were mentioned ear- 
lier by, e.g., Preston et al. (1965). Szeidl (1988) reported on 
weaker and stronger cycles in the star. From the Blazhko 
cycles covered by the data presented in this paper we also 
see variation, even from one Blazhko cycle to the next. Such 
a variation may be due to additional cycles in the star that 
are yet to be unraveled, or to non-periodic changes. The fact 
that we see residual signal around the position of the main 
frequency and its harmonics after subtracting the average 
frequency and the multiplet components is most likely due 
to these variations (periodic or transient). Finally, the large 
corrections needed to adjust the flux measurements during 
some phases (Fig. 3), particularly in Q2, are probably not 
critical with regard to the detected additional variations. As 
pointed out in Section [2~2| these corrections do not introduce 
significant "scatter" in the frequency results. Hence we can 
safely claim that most of the "scatter" can be attributed to 
variations in the Blazhko effect from cycle to cycle. 

RR Lyr itself is reported to have a longer cycle, of about 
four years, as first mentioned by Detre & Szeidl (1973). At 
the end of such a four-year cycle, the star's Blazhko effect 
is reported to weaken in strength, after which a new cycle 
starts with increased strength and with a shift (of about 
0.25) in the Blazhko phase tp. 

Kolenberg et al. (2006, 2008) list ephemerides for the 
Blazhko phase based on data gathered over the past decade. 
However, with a changing Blazhko period and large gaps in 
the data, it is nearly impossible to determine whether there 



has been a phase shift in the Blazhko effect at a distinct mo- 
ment in time. If the abrupt phase shift is a real phenomenon, 
we hope to be able to observe it in the Kepler data to come. 
Due to the period doubling it is not trivial to establish pre- 
cise ephemerides for RR Lyr, especially with long-cadence 
data. Forthcoming data gathered in short cadence will allow 
us to do so. 



4.4 Light curve changes and bump behaviour 

Fig.[9] also shows the periodic behaviour of the so-called 
"bump" in the light curve that appears just before mini- 
mum light. The bump occurs at an earlier phase around 
Blazhko minimum and later around Blazhko maximum. Also 
its strength shows a dependence on the Blazhko phase: it is 
most distinct during Blazhko minimum and becomes weaker 
around Blazhko maximum. A similar behaviour of the bump 
was also observed by Kolenberg et al. (2009) in the Blazhko 
star SS For. Guggenberger & Kolenberg (2007) investigated 
it more quantitatively for both SS For and RR Lyr. 

The bump has been explained by a collision between 
layers in the deep atmosphere, leading to shock waves 
that, depending on their direction, are explained by the 
"infall model" or the "echo model" (Hill 1972; Gillet & 
Crowe 1988). In their seminal work on spectroscopic data of 
RR Lyr, Preston et al. (1965) concluded that a displacement 
of the shock forming region occurs over the Blazhko cycle. 
This induces an earlier or later occurrence of the bump in 
the light curves. Our observations suggest that this collision 
may be stronger and earlier when the star is at minimum 
Blazhko phase. According to the findings by Jurcsik et al. 
(2009) for MW Lyr, this is also the phase when the star's 
mean radius is smaller (see their fig. 14), which could explain 
the earlier occurrence of the collision. Although the radiative 
hydrodynamical models of Stothers (2006) are very crude, 
especially regarding the amplitudes of small features, the 
phase of the bump in his models also occurs as observed 
when the pulsation amplitude is low. So it may be that this 
phasing is a consequence of an overall lower light amplitude, 
with the "echo" arriving at the surface sooner because it is 
running through a more compact star. 

Moreover, the so-called "hump" , the shoulder on the ris- 
ing branch of the light curve, only occurs at certain phases in 
the Blazhko cycle. In another early paper, Struve (1948) re- 
ported striking hydrogen emission line profiles in RR Lyr at 
particular, always "rising Blazhko" phases (increasing light 
amplitude of the cycle light curves). This happens to be 
the Blazhko phase interval when we see a more pronounced 
hump on the rising branch. The sudden slope change of the 
rising branch may be connected with layers colliding (hence 
emission in the spectra) and thus slowing down the motions 
(hence slope change in the light curve). With long cadence 
data the rising branch is covered with too few data points 
(often only three). Both the bump and hump characteris- 
tics can be investigated in detail with short-cadence data of 
RR Lyr. 

4.5 Constraints on the existing Blazhko models 

How do the Kepler observations of RR Lyr provide con- 
straints for improving the existing Blazhko models? The 
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most striking features we found, for the first time, in this 
data set are the following: 

(i) We found additional frequencies: quintuplet compo- 
nents and, most importantly, the period doubling phe- 
nomenon. 

(ii) We detected extreme light curve variations measured 
with high accuracy, e.g., the hump and bump feature. 

(iii) We observed that there is no strict repetition from 
one Blazhko cycle to the next, and that the Blazhko period 
may also vary on this short time scale. 

In addition, the most important findings so far from long 
cadence data of all 29 Kepler RR Lyrae stars are (Benko et 
al. 2010): 

(i) About half of the stars, 14 out of 29 studies RR Lyrae 
stars in our Kepler sample, are modulated. 

(ii) Period doubling and additionally excited higher-order 
radial overtones detected in seven stars, all of them Blazhko 
stars. 

(iii) All 14 Blazhko stars show both period modulation 
and amplitude modulation, in differing degrees. 

With increasingly impressive data sets and subsequent 
new and detailed findings, it is clear that the models so 
far proposed for the Blazhko effect are lagging behind. Any 
model for the Blazhko effect has to take into account both 
amplitude and phase modulation, as pointed out by Benko 
et al. (2010). Moreover, explanations that strictly impose a 
single modulation period, (e.g., the star's rotation period) 
not allowing for variation of the modulation, are definitely 
no longer viable. The moving light curve features (hump, 
bump) can provide additional clues about the stellar dynam- 
ics varying over the Blazhko cycle. Moreover, any plausible 
model has to be able to account for the (transient) occur- 
rence of period doubling and possible excitation of higher 
order overtones. 

The recent findings by Szabo et al. (2010) in connection 
with the period doubling suggest a 9:2 resonance between 
the radial mode and the 9th-order radial overtone (a strange 
mode) to be responsible for the period doubling. 

In the scenario proposed by Stothers (2006) a varia- 
tion of the star's convection is the cause for the Blazhko 
modulation. As the star's mean stellar parameters (temper- 
ature, radius and luminosity) vary over the Blazhko cycle 
(see Jurcsik et al. 2009), one can reasonably assume that its 
turbulent/convective structure also changes. It is unclear, 
however, what could lie behind this change. We note that 
transient magnetic fields at the origin of the variation in 
convective turbulence, as postulated by Stothers (2006), are 
most probably undetectable with present-day instrumenta- 
tion, and thus their existence would be hard to prove. 

It is clear that it is time to revise or expand the exist- 
ing models for the Blazhko effect, as well as explore alter- 
native explanations, using the constraints described above. 
Perhaps it is also time to revisit previously discarded sce- 
narios, such as models based on radial mode resonances (see, 
e.g., Goranskij et al. 2009; Kovacs 2009), or explore radial 
hydrodynamical models in much more detail. After the mod- 
els by Borkowski (1980) and Moskalik (1986) were dismissed 
by lack of confirmation from hydrodynamical models, most 
attention has been on resonances involving nonradial modes. 



5 CONCLUSIONS AND OUTLOOK 

The first releases of the Kepler data of the prototypical star 
RR Lyr are revolutionary because: 

• These data prove that bright stars such as RR Lyr can 
be observed with Kepler, as the saturated flux is conserved 
to a high degree (less than 2.5 x 10~ 3 fractional error even at 
maximum flux) . Generally the philosophy is not to waste any 
pixels on a star of which one cannot capture all the useful 
flux. RR Lyr is an unusual case due to its high variability. 
Because there are measurements in which the star is well 
captured and there is a lot of saturated flux in columns 
captured in every measurement, we can make good estimates 
of missing flux and thus recover photometry. Such unusual 
cases were not considered earlier in the Kepler Mission. 

• From the three Blazhko cycles covered by the data set 
(Q1+Q2 data), we obtain f = 1.76416±0.00001 d" 1 for the 
average radial pulsation frequency and /b = 0.0256 ±0.0002 
d _1 for the Blazhko frequency. The Blazhko period is Pb = 
39.1 ± 0.3 d, which confirms the earlier findings that the 
star's modulation period has gradually become shorter over 
the past decades (Kolenberg et al. 2006) . 

• There are, however, also indications that the Blazhko 
period is variable on shorter time scales, because our data 
clearly indicate that there is no strict repetition from one 
Blazhko cycle to the next. 

• We clearly detect quintuplet components in the fre- 
quency spectrum of the Kepler RR Lyr data. It is the first 
time these are found in the star. Their behaviour is similar 
to that of the triplet components. 

• Moreover, we find clear evidence for additional frequen- 
cies in the data. The most striking new frequency pattern 
is given by the so-called half-integer frequencies, that oc- 
cur at 2 2 1 /o, with the highest peaks around 3/2/o, 5/2/o, 
l/2/o, etc. Their occurrence is connected to the observation 
of alternating higher and lower maxima throughout (cer- 
tain phase intervals of) the Blazhko cycle. This phenomenon, 
called period doubling, has been observed for the first time 
in Kepler data (Kolenberg et al. 2010a). Meanwhile, we have 
observed it in several Blazhko stars observed with Kepler, 
and its presence may be connected to the Blazhko effect. The 
period doubling is variable over the course of the Blazhko 
cycle, and also does not repeat identically from one Blazhko 
cycle to the next (Szabo et al. 2010). 

• The Fourier parameters derived from individual cycles 
oscillate around the Fourier parameters derived from several 
consecutive cycles, for which the period doubling effect is 
smoothed out. 

• The increase of the Fourier parameter R21 with decreas- 
ing pulsation amplitude Ai is counterintuitive, because pul- 
sations with lower amplitudes (lower Ai) are expected to be 
less nonlinear (lower i?2i). This observation provides con- 
straints for modelling the physics of the Blazhko modula- 
tion. 

• The position and strength of the bump and hump vary 
over the Blazhko cycle. A detailed study of these features 
in the pulsation cycle, in photometric as well as spectro- 
scopic data, may shed new light upon the pulsation and 
shock propagation in RR Lyrae stars and the understand- 
ing of the Blazhko effect. 
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After the Q2 observations, it was realized that RR Lyr's 
Kepler KIC magnitude was not bright enough to match 
the observations, and it was adjusted upwards. As a con- 
sequence, the standard aperture assignment algorithms as- 
signed it to an aperture that went off silicon. Thus the tar- 
get was rejected, and was not observed in Q3 and Q4. The 
first Kepler findings described in Kolenberg et al. (2010a), 
however, triggered further interest in the star. Therefore, a 
custom aperture was devised for RR Lyr (by Steve Bryson), 
reducing the amount of originally assigned pixels from 433 to 
about 150. Thanks to the calibration work described in this 
paper, RR Lyr will likely be scheduled to be a target of Ke- 
pler for the rest of the mission. At the time of writing of this 
paper, RR Lyr is being observed through the Kepler Guest 



Observer (GO) program (http://keplergo.arc.nasa.gov/) 



As further Kepler data of RR Lyr become available, we 
will be able to perform more focused analyses: 

• A longer time base can lead to a better resolved fre- 
quency spectrum. However, due to the period variation and 
the transient phenomena, the frequency spectrum will be 
smeared out. Hence, in reality we will probably not obtain 
"sharper" peaks with a longer time base, but several closely- 
spaced peaks. New and additional Kepler data will enable us 
to analyse the temporal variation of the frequencies. They 
will also help uncover whether the variation between consec- 
utive Blazhko cycles is due to additional cycles in the star, 
or to transient or quasi-periodic changes. 

• Kepler short-cadence data will allow us to study short- 
time effects and light curve features such as the hump and 
bump. Simultaneous ground-based data are planned for the 
period of the short-cadence Kepler observations. 

• With Kepler's potential coverage of RR Lyr's light vari- 
ations over 3.5-5 y, many mysteries still surrounding this 
prototypical star can be solved. We may witness the start 
of a new 4-year cycle (Detre & Szeidl 1973), as RR Lyr is 
known to have. 

• With Kepler we can anticipate a dramatic overhaul in 
the models for the Blazhko effect. The constraints provided 
by the ultra-precise data of RR Lyr motivate us to revisit, 
revise or expand the existing models for the Blazhko effect, 
as well as explore alternative explanations. 
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